Overcoming non-Markovian dephasing in single photon sources through post-selection 
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We study the effects of realistic dephasing environments on a pair of solid-state single-photon 
sources in the context of the Hong-Ou-Mandel dip. By means of solutions for the Markovian 
or exact non-Markovian dephasing dynamics of the sources, we show that the resulting loss of 
visibility depends crucially on the timing of photon detection events. Our results demonstrate that 
the effective visibility can be improved via temporal post-selection, and also that time-resolved 
interference can be a useful probe of the interaction between the emitter and its host environment. 
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Indistinguishable photons lie at the heart of linear opti- 
cal [l| and distributed light-matter quantum information 
processing (QIP) Q. In particular, interference effects 
such as photon bunching and 'which-path' erasure are 
central to a range of QIP protocols, leading to a recent 
drive for highly efficient single photon sources. Solid- 
state systems such as quantum dots (QDs), individual 
molecules, or colour centres in diamond [3| are emerging 
at the forefront of this new technology due to the po- 
tential they offer for realizing not only efficient but also 
controllable, deterministic single photon devices. 

However, solid-state systems are inherently noisy and 
thus their optical transitions are typically subject to en- 
ergy fluctuations, for example due to coupling to phonon 
modes in the host system. These fluctuations tend to re- 
duce the phase coherence of the emitted light, and hence 
degrade the visibility of interference effects [J] . Dephas- 
ing can be viewed as broadening of the optical transition, 
which reduces the indistinguishability of the emitted pho- 
tons. In the case of QIP schemes 0, 0], this dephasing 
leads to logical errors in the calculation. Overcoming 
this effect is important in view of the practical applica- 
tions of single photon sources. The realization of Fourier 
transform (i.e. lifetime) limited single photon devices 
has therefore become the goal for an extremely active 
area of research p|. Moreover, understanding the inter- 
action between the emitter and its host environment is 
an interesting problem in its own right. 

In this Rapid Communication we propose that these 
issues can be addressed with the aid of temporal post- 
selection of the photons emitted from a noisy source. If 
sufficiently fast photodetectors are available, it is possible 
to post-select a sub-ensemble of detection events based on 
their detection times. By considering a Hong-Ou-Mandel 
(HOM) type two-photon interference setup (see inset of 
Fig. [T|), we show that such post-selection can substan- 
tially improve the effective visibility in an interference 
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experiment, albeit at the expense of effective source effi- 
ciency. This allows significant improvement in the fidelity 
of QIP protocols, at the expense of success probability, 
which nevertheless will be of benefit in near- future exper- 
iments with noisy sources. Furthermore, we show that 
post-selection can yield important information about the 
basic physics of the interaction between the emitter and 
its host environment. By measuring the HOM visibility 
as a function of time, one can directly probe the dynamics 
of the system-environment interaction. 

We make use of an explicit model of a single photon 
source which is coupled both to the modes of the optical 
field, and to a solid-state environment. While the optical 
modes can be observed via single photon detectors, the 
experimentalist has no direct access to the environment 
degrees of freedom. This is a familiar situation in quan- 
tum optics, which could be dealt with using the standard 
technique of quantum trajectories [5|, whereby both the 
system-photon and system-environment couplings can be 
treated pertubatively within the Born-Markov (BM) ap- 
proximation 0, H, However, although the BM ap- 
proximation is typically very well satisfied for the system- 
photon coupling, in general, this approximation may 
break down for the system-environment coupling, and 
non-Markovian effects can be significant [![. Here, we 
move beyond the BM approximation for the solid-state 
environment, and find an exact analytic solution for the 
dynamics induced by the system-environment coupling. 
We treat the system-photon coupling within the BM ap- 
proximation, which allows us to compute the observed 
photon statistics in a straightforward manner. 

We model our single photon sources as a pair of two- 
level systems each with ground state \g), excited state 
|e), and equal energy splitting E. Interaction with the 
radiation field induces spontaneous decay from state |e) 
to \g) with associated photon emission. In addition, both 
systems couple to their own bath of harmonic oscillators, 
representing the solid-state environment. We assume the 
sources are placed far enough apart that no interactions 
occur between them, either direct or mediated by the 
solid-state baths. A relevant physical implementation 
could, for example, be two optically active self-assembled 
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QDs @ or nitrogen- vacancy (NV) centres [l(| , each cou- 
pled to their own cavity and operating within the Pur- 
cell regime 0, EH- Emitted photons are interfered on 
a 50:50 beam splitter, without any delay introduced be- 
tween the arrival times, and subsequent detection occurs 
at the photon counters D+ and D-. We will not treat im- 
perfections in the setup such as limited photon detection 
and collection efficiency, or system frequency mismatch, 
in order to concentrate on the impact of the dephasing 
baths. However, in principle, such effects could be in- 
corporated into the formalism fl|. In fact, the only dif- 



ference due to realistic detection efficiency would be a 
reduction in the total number of two-photon events; if 
one were to post-select only the sub-ensemble of events 
where two photons were observed from each excitation, 
one would obtain exactly the results presented here. 

The Hamiltonian describing our setup is given by 
H = H Q + Hi, with H = £\ H jfi = \e) j (e\ + 

Hj,B + Hj^x) representing each two-level system and 
the harmonic oscillator baths in isolation (j — 1,2 de- 
notes the two qubits). Here, Hjb — Ek w 3> a J k a J, k 

and Hj t x = E q Vj.qbj q &j, q describe the solid-state and 
photon environments, with frequencies ujj^ and ^-, q re- 
spectively (h = 1 throughout). The interaction term is 
given by Hj = T,j H j,i = Ej(-%SB + #i,sx), where 

Hj,SB = < , < zLk^Z/ k"; k • Wj.k'O.k* and H j,SX = 

E q (lff)j (e\ /i,q 6 ],q + _ \ e )j (fl /j, q & i,q) account, respec- 
tively, for system-solid-state bath interactions via the 
excited states |e) ■ with coupling gj^, and for system- 
photon interactions (within the dipole and rotating-wave 
approximations [13]) through the raising and lowering 
operators |e)^ (g\ and \g)j (e|, with coupling / Jiq . The 

bath creation (annihilation) operators at k (a^k) and &t q 
(bj. M ) obey boson statistics, and none of the baths share 
any common modes. We note that the form of system- 
solid-state bath interaction (-ff^se) assumed above is par- 
ticularly relevant to the study of short-time dephasing 
behaviour [l3| , such as that observed experimentally due 
to phonons in QDs 

Moving into the interaction picture with respect to Ho, 
the exact evolution of the combined system-environment 



density operator xi(^) 



*i = -i[Hi(t),x(0)] 



ot x(t)e 



-iH t 



is given by [U 



dT[Hl(t),[Hl(T),Xl(T)}}, (I) 



where H(t) = e^Hje-*"** = Ej^'.sbW + ^sxW). 
We are interested here in the regime in which the system- 
solid-state bath dynamics is slow on the timescale set by 
the radiation field correlation time (which is typically 
extremely short, tv. ~ 1/E ~ 1 fs for E = 1 eV [l3|). 
In this case, we may perform a Born approximation 
on the density operator factorizing it into a system- 
solid-state bath contribution p and a radiation field part 
R as x = p ® Ri ® R2, assuming the baths are ini- 
tially uncorrelated. It is important to stress that we do 
not further factorize p and so fully account for correla- 



tions built up between the system and solid-state bath 
during their combined evolution. Inserting this form 
for x m f° Eq. and tracing over the photon baths 
we find the simpler form pi = —iJ2j[Hj-SB(t),Pi(t)] — 

Ej Jo dTTl jM H j,sx(t),[ H j,sx.( T ),Pi( T ) R j]]: where we 
assume that Rj are thermal equilibrium states. The in- 
terpretation of this equation is clear; the first term de- 
scribes the exact evolution of pi due solely to iJi,sB an d 
-H2,sb, while, within a Markov approximation, the sec- 
ond term gives rise to the usual Lindblad operators asso- 
ciated with spontaneous emission and defines the transi- 
tion rates jj for each qubit [!, [Hj] ■ 

The quantum jump formalism is therefore particularly 
well suited to studying the dynamics of our system, here 
subject to continual observation by detectors D±. We 
define jump operators corresponding to photon emission 
and detection as Cj — yfyj\g)j (e|, projecting each two- 
level system from its excited to ground-state (for sim- 
plicity, we take 71 = 72 = 7)- The action of the beam 
splitter is to mix these operators, erasing the informa- 
tion on where the emis sion originated, and leading to 

c± = (ci ± c 2 )/V2 = v / r72(lfl')i ( e l ± 1^)2 ( e l); as J um P 
operators for a photon count in D± , respectively. Assum- 
ing the initial state of the sources to be |e) x (8> |e) 2 = |ee), 
and that all emitted photons are collected and detected, 
within any increment of time dt the system may evolve 
in one of only two possible ways. If no photon is ob- 
served, we know that neither source has returned to its 
ground-state. In this case, the combined system-solid- 
state bath evolution is generated by the conditional no- 
jump master equation pi = — i(H e g(t)pi — piH^ s (t)), with 
non-Hermitian effective Hamiltonian 

H cS (t) = Y, H iMt) - {i/2)c\c + - (i/2)c f _c_. (2) 
j 

We define the corresponding non-unitary, no-jump time- 
evolution operator by 



Ui(t h tf) = T exp (~*J t 1 dtH eS (t) 



(3) 



where T is the time-ordering operator. The state gen- 
erated by jh(tf) = Ui(ti,tf)pj(ti)Uj(ti,tf) will be un- 
normalised and Tr(/>i(i/)) can thus be interpreted as the 
probability that no detector clicks are observed within 
the time interval to tj. Accounting for the effect of the 
time-ordering operator in Ui(U,tf) is straightforward if 
we employ a Magnus expansion of the exponential [l5| , 
since all terms of the form [H e g(t), [H e e(t'), H e g(t")]] and 
higher disappear, allowing us to solve for the exact dy- 
namics due to the solid-state environment [l6| . 

A click in either detector, on the other hand, implies a 
discontinuous lowering of the system state and we apply 
the appropriate jump operator 



pi(t + dt) = c±pi(t)c'±/Tr(&±c±pi(t)), 



(4) 



for a count in D±. Here, Tr(cJ_c±pi(t)) is the probability 
that a click occurs in detector D± in the interval t to 
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FIG. 1: Main: Visibility against lu c t for: an ohmic bath 
with r„ =1 (T) = Aln(l + w 2 -r 2 ) + A In (sinh(7rr//?)/(7rr//3)) 
(blue, solid line); a superohmic bath with Y n =s{T) = 
A^r 2 (3 + o, 2 r 2 )/(l + lu 2 t 2 ) 2 + A(tt 2 /3(uj c f3) 2 - 1/(lj c t) 2 + 
ir 2 csch 2 (tvt/ (3)/ (i0cP) 2 ) (yellow, dashed line); and Markovian 
dephasing of rate Tm(t) = Anr//3 (red, dotted line). Param- 
eters: A = 0.5, u>cP = 10. Inset: Schematic of the interference 
setup. Two photons incident on a 50:50 beam splitter (BS) 
interfere such that they are always detected at the same out- 
put, provided the photons themselves and the paths they take 
are identical. The resulting absence of coincidence counts be- 
tween the detectors is known as the HOM dip. Its visibility 
measures the indistinguishability of the incident photons. 



t + dt given that that the system is in a normalized state 
pi(t) at time t. 

Eqs. (J2EJ) allow us to calculate the full time-resolved 
statistics of the photodetector counts, and thereby as- 
sess the impact of source dephasing on the visibility of 
the HOM dip. The sources are again assumed to be ini- 
tialised in state |ee), with the baths in thermal equilib- 
rium, p jjB = e~ /3H ^ B /Tr(e~ /3ff ^ B ), where (3 = l/k B T, k B 
being the Boltzmann constant. The combined system- 
bath density operator is thus of initial form p(0) 
1 1 (ee| <g> pi t B <8> P2,b- It is instructive to first calcu- 
late the intermediate state of the system, ps, at time i 2 , 
conditional on observing the first detector click at time 
ti < ti in detector D + , but before the second click. From 
Eqs. on tracing over both environments, we find 



~7 T r 



PS 



\ge) (ge\ + \eg) (eg\+ e -^ +r ^ 



x (e***^) \ge) (eg\ + e -*(*i.*»> \eg) (ge\)] , (5) 



where r = i 2 — t\ is the time elapsed since the first 
click. Here, T 3 {t) = £ k \a 3 , k (*i, t)/V2\ 2 coth (p Uj , k /2), 
with a j:k (t l7 T) = (5 i:k /w Jik )e^- ktl (l-e l ^- kT ), are time- 
dependent decoherence functions, characterising the un- 
desirable effects of the baths on the system evolution. 
The phase factor is 0(ii,i 2 ) = A 2 (ti,t 2 ) — Ai(ii,i 2 ), 

where Aj-(ii,i 2 ) = \9jM/ UJ j,k\ 2 ( UJ jM T + 2sinwj',k*i - 
2sinujj^t2 + sina>j ;k r). We restrict ourselves to the case 
when both sources couple equally to the same form of en- 
vironment, giving Ti(t) — r 2 (r) = T(t) and (/>(ii,i 2 ) 
I). This is realized experimentally when two consecutive 
photons are interfered from the same source Q. 

The probability density of the second click being ob- 
served in either the same detector D + (= p(i 2 , +\ti, +)) 



or opposite £>_ (= p(i 2 , — \t±, +)) during a time inter- 
val i 2 to i 2 + dt, conditional on the first click occur- 
ring in D+ at time t%, may now be calculated from 
p(t 2 ,±|*i,+) = Tr(4c±p s (r)). We find 



p(t 2 ,±\h,+)=p(±\ 



7T 



lie" 



-ar(r) 



(6) 



which reduces to t) = r ye~' yT and p(— 1+, r) = in 

the absence of dephasing, then implying that both pho- 
tons exit at the same port of the beamsplitter (i.e. perfect 
HOM interference). Dephasing of the sources therefore 
destroys the interference effect, and by implication the in- 
distinguishability of the emitted photons, by an amount 
that depends not only on the form and strength of the 
system-bath interactions, but also on the difference in 
detection times r. This may be seen more clearly by con- 
sidering the visibility of the HOM dip, which we define 
in the ideal detector limit to be 



y(r) 



P (+\+,t)-p(-\+,t) 



p(+\+,t)+ P (-\+,t) 



-2r(r) 



(7) 



measuring the normalised difference in conditional detec- 
tion probabilities for D + and D_ at time t 2 . In the ab- 
sence of any dephasing v(r) = 1 for all r as we consider 
no further imperfections in the setup, while in general 
v(t) is completely determined by the decoherence func- 
tion r(r). There is no dependence on the spontaneous 
emission rate at the perfectly time-resolved level. 

We plot the visibility as a function of detector click 
separation in Fig. [1] To illustrate the general behaviour, 
we cast the decoherence function in terms of the spectral 
density [H, J(u) = £ k \g k \ 2 S{u; - cj k ) = Auj n e^^% 
as r(r) = J™(J(w)/uj 2 )(1 - cos wt) coth (/3w/2). Here, 
ui c is a high frequency cut-off and A an appropriate cou- 
pling strength. We consider both ohmic (n — 1) and 
superohmic (n = 3) spectral densities as they are rele- 
vant, for example, to carrier-phonon interactions in the 
solid-state 16]. We also plot the performance with a 
Markovian rate for comparison, where Tm(t) is given 
by the thermal limit of the ohmic case. We see that 
for both the ohmic and Markovian cases the visibility 
tends asymptotically to zero as the click separation be- 
comes large, although the ohmic behaviour differs qual- 
itatively on shorter timescales due to non-Markovian ef- 
fects. In the case of a superohmic spectral density the 
evolution is entirely non-Markovian, with remnant visi- 
bility ^ n=3 (oo) = exp (-2A(1 + 7r 2 /3(w c /3) 2 )). 

Clearly, for ideal detectors of time resolution much 
shorter than any system-bath dynamics, Fig. [T] sug- 
gests that by post-selecting detection events closely 
spaced in time [181 ] the visibility of the two-photon in- 
terference may be improved, overcoming noise inher- 
ent within the sources. In fact, this is also true for 
non-identical sources, where the visibility becomes v = 
e -(ri(r)+r 2 (r))| cos ^ i)ti +r )|_ Her6j post-selection 

will still improve the quality of the interference, though 
now it is not simply r that determines the loss in visi- 
bility, but also the detection time t-y itself. Although the 



4 




5 10 15 20 

o) c A 



FIG. 2: Visibility against ui c A for the ohmic (blue, solid line), 
superohmic (yellow, dashed line), and Markovian (red, dotted 
line) baths. Parameters are as in Fig. [T] with 7/cj c = 0.01. 

new factor in the visibility is oscillatory, in practice it is 
impossible to measure the detection time with arbitrar- 
ily high precision. Hence, we propose setting a detection 
time post-selection window. In this situation, for non- 
identical sources, it is generally advantageous not only 
to post-select closely spaced detection events, but those 
in which the detection time t\ is short as well. 

Returning to the case of identical sources, we illus- 
trate this idea in Fig. [2j where a modified visibility func- 
tion ✓(A) = |(p++(A) -p+_(A))/(p++(A) +p+_(A))| 
is plotted. Here, p ++ (A) — f Q c drp(+\+,T) and 

p^ (A) = Jq" A drp(— |+, r) are the relevant conditional 

probability densities now integrated over the window 
width oj c A. We see that, in principle, the improvement 
due to post-selection can be extremely good since the vis- 
ibility will always approach unity as A — ► 0. However, 
any improvement necessarily comes at the cost of effec- 
tive source efficiency and practical considerations on the 
number of photons that can reasonably be discarded lim- 
its what would be possible for any particular setup. It is 
important to stress that such an improvement can only be 
observed in a time-resolved experiment [l8[ , though even 
if this resolution is fairly poor some improvement should 
be possible when there is a long-time decay of visibility 
towards a fixed value, as in the ohmic and Markovian 
examples considered here. In contrast, in the 'bad detec- 
tor' limit of zero resolution, we must integrate the condi- 
tional probability densities over all r. Here, assuming a 
set system-bath interaction strength, the simplest way to 
improve the visibility is to engineer a larger effective de- 
cay rate such that the system emits on timescales shorter 
than those of the dominant dephasing processes [l9| . 



In the opposite limit, when the system spontaneous 
emission rate is small compared to the typical rate for 
dephasing Q, the form of the system-bath interaction 
can itself be probed with fast detectors by successively 
widening the post-selection window (scanning the visibil- 
ity in time as in Fig IS}. Provided that the post-selection 
window size is kept much shorter than the excited state 
lifetime (7 A < 1), then the visibility has essentially no 
dependence on the value of 7 itself, allowing the bath- 
induced dynamics to be isolated and characterised in- 
dependently from spontaneous decay processes. For ex- 
ample, in this regime the visibility is well described by 
v' M (A) = /?(l-e" 2A7rA//3 )/2AvrA in the Markovian case, 
and by v' n=1 {A) = B_ A a (1/2,1 - 2A)/2iA in the low 
temperature ohmic case, where B z (a, b) is the incomplete 
Beta function. Such measurements therefore offer the po- 
tential for clear signatures of deviations from Markovian 
dephasing behaviour, as well as a means to assess system- 
bath coupling strengths. 

To summarise, we have shown analytically how to in- 
corporate both Markovian and non-Markovian dephasing 
baths into the quantum jump formalism for emitting two- 
level systems. We have analysed the impact of such baths 
on HOM interference from a pair of solid-state single pho- 
ton sources, demonstrating that the dephasing may be 
overcome by post-selection providing that sufficient time- 
resolution exists within the photodetectors. In general, 
time-resolved two-photon interference can provide signif- 
icant insight into solid-state decohcrence processes. As 
fast detectors down to sub-picosecond time resolution are 
now available [20| , much shorter than the typical excited 
state lifetimes of nanoseconds or greater in, for example, 
QDs or NV centres, experiments to explore such phe- 
nomena are feasible with current technology in a number 
of different physical systems, and will become more so 
as technology improves. Furthermore, it should also be 
possible to apply our formalism to gauge the impact of 
dephasing on a wide range of distributed entangling pro- 
tocols that exploit photon interference effects. 
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